Comparison of two quantum-cluster approximations 
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We provide microscopic diagrammatic derivations of the the Molecular Coherent Potential Ap- 
proximation (MCA) and Dynamical Cluster Approximation (DCA) and show that both are $- 
derivable. The MCA (DCA) maps the lattice onto a self-consistently embedded cluster with open 
(periodic) boundary conditions, and therefore violates (preserves) the translational symmetry of the 
original lattice. As a consequence of the boundary conditions, the MCA (DCA) converges slowly 
(quickly) with corrections 0(1/L C ) (0(1/L^.)), where L c is the linear size of the cluster. These 
analytical results are demonstrated numerically for the one-dimensional symmetric Falicov-Kimball 
model. 
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Introduction One of the most active areas in con- 
densed matter physics is the search for new methods to 
treat disordered and correlated systems. In these sys- 
tems, especially in three dimensions or higher, approxi- 
mations which neglect long ranged correlations are gener- 
ally thought to provide a reasonable first approximation 
for many properties. 

Perhaps the most successful of these methods are the 
Coherent Potential Approximation (CPA) jlj and the Dy- 
namical Mean Field Approximation (DMFA)^], ||, ||, 
||, for disordered and correlated systems, respectively. 
Although these approximations have different origins, 
they share a common microscopic definition. Both the 
DMFA(| and the CPA@ may be defined as theories 
which completely neglect momentum conservation at all 
internal diagrammatic vertices. When this principle is 
applied, the diagrammatic expansion for the irreducible 
quantities in each approximation collapse onto that of a 
self-consistently embedded impurity problem. 

Many researchers have actively searched for a technique 
to restore non-local corrections to these approaches. 
Here, we discuss just two approaches which are fully 
causal and self-consistent: the Molecular Coherent Po- 
tential Approximation (MCA)fl 13] and the Dynamical 
Cluster Approximation (DCA)|j, [H)| [n], 0. Recently 
the Cellular Dynamical Mean Field Approach || was pro- 
posed for ordered correlated systems, while the Molecular 
Coherent Potential Approximation has traditionally been 
applied to disordered systems. Since both methods share 
a common microscopic definition we use the term MCA 
to refer to both techniques in the following. 

While the MCA is traditionally defined in the real space 
of the lattice, the DCA is traditionally defined in its re- 
ciprocal space. In the MCA, the system lattice is split 
into a series of identical molecules. Interactions between 
the molecules are treated in a mean-field approximation, 
while interactions within the molecule are explicitly ac- 
counted for. In the DCA, the reciprocal space of the 
lattice is split into cells, and momentum conservation is 
neglected for momentum transfers within each cell while 
it is (partially) conserved for transfers between the cells. 
These approximations share many features in common: 
they both map the lattice problem onto that of a self- 
consistently embedded cluster problem. Both recover the 
single site approximation (CPA or DMFA) when the clus- 



ter size reduces to one and become exact as the cluster 
size diverges. Both are fully causal 0, Q, and provided 
that the clusters are chosen correctly [pT they maintain 
the point group symmetry of the original lattice problem. 
Here, we provide a microscopic diagrammatic derivation 
of both the MCA and the DCA, and explore their con- 
vergence with increasing cluster size. 

Formalism For simple Hubbard- like models, momen- 
tum conservation at each vertex is completely described 
by the Laue function 



,ix.(ki+k2-k3-k 4 ) _ jyr 



-k 2 ,k 3 +k 4 



(1) 



where ki, k2 (k3, ki) are the momenta entering (leav- 
ing) the vertex. Muller-Hartmann || showed that the 
Dynamical Mean Field (DMF) theory may be derived by 
completely ignoring momentum conservation at each in- 
ternal vertex by setting A = 1. Then, one may freely 
sum over all of the internal momentum labels, and the 
graphs for the generating functional $ and its irreducible 
derivatives, contain only local propagators. 

The DCA and MCA techniques may also be defined by 
their respective Laue functions. Since our object is to 
define cluster methods we divide the original lattice of N 
sites into N/N c clusters (molecules), each composed of 
N c = sites, where D is the dimensionality. We use 
the coordinate x to label the origin of the clusters and 
X to label the N c sites within a cluster, so that the site 
indices of the original lattice x = X + x. The points x 
form a lattice with a reciprocal space labeled by k. The 
reciprocal space corresponding to the sites X within a 
cluster shall be labeled K, with K a — n a ■ 2ir/L c and 
integer n a . Then k = K + k. Note that e* K x = 1 since 
a component of x must take the form m a L c with integer 



In the MCA, we approximate the Laue function by 



y e iX-(K 1 +K 2 -K 3 -K4-|-ki-|-k 3 -k3-k4) 



(2) 



Thus, the MCA omits the phase factors e lkx resulting 
from the position of the cluster in the original lattice 
but retains the (far less important) phase factors e lk x 
associated with the position within a cluster. In the DCA 
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we also omit the phase factors e lkx , so that 

Adc — ^c < 5ki+K2,K3+k 4 • (3) 

Both the MCA and DCA Laue functions recover the ex- 
act result when N c — > oo and the DMFA result, A = 1, 
when N c = 1. 

If we apply the MCA Laue function Eq. || to diagrams 
in <f>, then each Green function leg is replaced by the 
MCA coarse-grained Green function 

G(Xi,X 2 ;x = 0) = 

e i ( K 1 +ki)-x 1G(Kl)K2 .k lj k 2 )e-^+^)- x ^ = 



N 2 ■ 



Kl,K 2 

ki,k 2 



N 2 



G(X 1 ,X 2 ,k 1 ,k 2 ) 



(4) 



ki,k 2 



or in matrix notation for the cluster sites Xi and X 2 
N, 



(5) 



The summations of the cluster sites X remain to be per- 
formed. Note that the inclusion of the phase factors e zkx 
in the MCA Laue-function Eq. |^ leads directly to a clus- 
ter approach formulated in real space that violates trans- 
lational invariance. Therefore the Green function is a 
function of two cluster momenta Ki, K 2 or two sites Xi, 
X2 respectively. 

If we apply the DCA Laue function Eq. ||, Green func- 
tion legs in $ are replaced by the DCA coarse grained 
Green function 



G(K) = §^G(K,k), 



N 



(6) 



since Green functions can be freely summed over the k 
vectors within a cell about the cluster momentum K. 
(We have dropped the frequency dependence for nota- 
tional convenience.) As a result, <& is a functional of the 
coarse grained Green function G(K) and thus depends 
on the cluster momenta K only. 

To establish a connection between the cluster and the 
lattice we minimize the lattice free energy 



F = -k B T(<f> c -tr [EG] +trln [G]) 



(7) 



where $ c is the generating functional calculated with the 
coarse-grained propagators, E is the lattice self-energy 
and G is the full lattice Green function. The trace in- 
dicates summation over frequency, momentum and spin. 
As we have discussed elsewhere, only the compact part 
of the free energy, <I>, is coarse-grained. F is stationary 
with respect to G when ^ = 0. This happens for the 
MCA if we estimate the lattice self energy as 

E(Ki,K a ;ki,k2) = 

]T e - l ( Ki+ki > Xi s MC (x 1 ,x 2 ) e l ( K2+k2 )- X2 . (8) 

Xi ,x 2 



Thus, the corresponding lattice single-particle propaga- 
tor reads in matrix notation 



G(k,z)= zI-e(k)-E MC (z) 



(9) 



where the dispersion e(k) and self-energy Emc(z) are 
matrices in cluster real space with 



[e(k)] Xl x 2 - c(Xi-X 2 ,k) 
1 

" N. 



(10) 



^ e i(K+k)(X! 
K 



-X, 



c K+k 



being the intracluster Fourier transform of the disper- 
sion. For the DCA, S(k) = Ebc'(K) is the proper ap- 
proximation for the lattice self energy corresponding to 
<&dc- The corresponding lattice single-particle propaga- 
tor is then given by 



G(K,k; z) 



r K+k 



Edc(K,z) 



(11) 



Both the MCA and DCA are optimized when we equate 
the lattice and cluster self energies. A similar rela- 
tion holds for two-particle quantities. Thus, with few 
exceptions jl3j, only the irreducible quantities on the clus- 
ter and lattice correspond one-to-one. 

The MCA (DCA) algorithm, with steps A^D, follows 
directly: A we first make an initial guess for the cluster 
self-en erg y matrix E. B This is used with Eqs. [| and ^ 
(^ and [ll]) to calculate the coarse-grained Green function 

G. C The cluster excluded Green function Q = [G + 
Emc]" 1 (0(K) = [G(K)- 1 + ^(K)]- 1 ) is defined to 
avoid overcounting self energy corrections on the cluster. 
It is used to D compute a new estimate for the cluster 
self-energy which is used to reinitialize the process. Once 
convergence is reached, the irreducible quantities on the 
cluster may be used to calculate the corresponding lattice 
quantities. 

In order to compare the character of the two different 
cluster approaches as a function of the cluster size N c it 
is instructive to rewrite the corresponding coarse grained 
Green-functions Eqs. [s]and ||to suitable forms by making 
use of the independence of the self-energy E on the inte- 
gration variable k. For the MCA coarse grained Green 
function we find 



G(z)= zl - e -T, M c(z) -Tmc(z) 



(12) 



with the "cluster-local" energy e = N c /N For 
the DCA we obtain a similar expression 



G(K, z ) = [z-e K - Edc(K, z) - T DC (K, z)\ 



(13) 



with the coarse grained average £k = N c /N^^. e(K, k). 
The hybridization functions T MC / DC (z) describe the 
coupling of the cluster to the mean field representing the 
remainder of the system. 
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The behavior of T for large N c is important. For the 
MCA, r averaged over the cluster and frequency 



o 



Xl ,X2 



2D 



(14) 



where L c = N c is the linear cluster size. A detailed 
derivation of this form will be published elsewhere. How- 
ever, since in the MCA the cluster is defined in real space 
with open boundary conditions, this form is evident since 
only the sites on the surface oc 2D ■ L®^ 1 of the cluster- 
couple to the effective medium. For the DCA we have 
previously shown that T(K) - 0(1/nI /D ) Jl| so that 
we obtain for the average hybridization of the DCA clus- 
ter to the effective medium 



K 



dc(K) ~ O 



(15) 



The DCA coarse graining results in a cluster in K-space; 
thus, the corresponding real space cluster has periodic 
boundary conditions, and each site in the cluster has the 
same hybridization strength T with the host. 

In both the MCA and the DCA, the average hybridiza- 
tion strength acts as the small parameter. The approxi- 
mation performed by the MCA (DCA) is to replace the 
lattice Green function G(k) = [zl — e(k) — S(k, z)]^ 1 
(G(K, k, z) = [z — e K+ k — X(K, k, z)] -1 ) by its coarse 

grained quantity G (G(K)) in diagrams for the self- 
energy S. Once the sums over k are performed, all terms 
which are lower order in 1/L C than T vanish. Thus the 
MCA (DCA) is an approximation with corrections of or- 
der f ~ 0(1/L C ) (~ 0{l/LD). 

Numerical Results To illustrate the differences in con- 
vergence with cluster size iV c we performed MCA and 
DCA simulations for the symmetric one-dimensional 
(ID) Falicov-Kimball model (FKM). At half filling the 
FKM Hamiltonian reads 

H = -tJ2(4d i+1 + h.c.) + Uj2(nf - l/2)(»/ - 1/2) , 



(16) 

with the number operators nf = d\d i and n{ = fjf i and 
the Coulomb repulsion U between d and / electrons re- 
siding on the same site. The FKM can be considered as a 
simplified Hubbard model with only one spin-species (d) 
being allowed to hop. However it still shows a complex 
phase diagram including a Mott transition and Ising-like 
charge ordering with the corresponding transition tem- 
perature T c being zero in ID. The dispersion in (ID) 
€k = 2icosfc; thus for t = 1/4 the bandwidth W — 1 
which we use as unit of energy. To simulate the effective 
cluster models of the MCA an the DCA we use a quan- 
tum Monte Carlo (QMC) approach described inQ. 

To check the scaling relations Eqs. [l4| and 15 
m Fig.|l| the average hybridization functions T 



we show 
mc and 

Trie f° r t ne MCA and DCA respectively at the inverse 
temperature (3 — 17 for U — W = 1. For N c = 1 both 



approaches are equivalent to the DMFA and thus f mc = 
Tdc- For increasing N c Tmc can be fitted by 0.3361/^ 
and f DC by 1.1946/A^ when iV c > 2. Cluster quantities, 
such as the self energy and cluster susceptibilities, are 
expected to converge with increasing N c like f . This is 
illustrated in the inset for the staggered (Q = it) charge 
susceptibility Xc(Q) of the cluster. 
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FIG. 1: The average integrated hybridization strengths T of 
the MCA (squares) and DCA (circles) versus the cluster size 
iV c when /3 = 17 and U = W = 1. The solid and dashed 
lines represent the fits 1.1946/iV,? and 0.3361/iV c respectively. 
Inset: Convergence of the cluster charge susceptibility for Q — 
7r. The solid and dashed lines are quadratic and linear fits, 
respectively. 

Since only the compact parts represented by <& of the 
lattice free energy (Eq. 0) are coarse-grained, this scaling 
is expected to break down when lattice quantities, such 
as the lattice charge susceptibility, are calculated. The 
susceptibility of the cluster Xc(Q) cannot diverge for any 
finite N c ; whereas the lattice x(Q) diverges at the tran- 
sition temperature T c to the charge ordered phase. Note 
that the residual mean field character of both methods 
can result in finite transition temperatures T c > for 
finite N c < oo. However as N c increases, this resid- 
ual mean field character decreases gradually and thus 
increased fluctuations should drive the solution to the 
exact result T c = 0. 

In the DCA|ll|], x(Q) is calculated by first extract- 
ing the corresponding vertex function from the cluster 
simulation. This is then used in a Bethe-Salpeter equa- 
tion to calculate x(Q)- T c is calculated by extrapolating 
xiQ) 1 to zero using the function xiQ) 1 (T — T c ) 7 
(see inset to Fig.||). This procedure is difficult, if not im- 
possible, in the MCA due to the lack of translational 
invariance. Here, we calculate the the order parame- 
ter m(T) = 1/N C m the symmetry broken 
phase. T c is then obtained from extrapolating m(T) to 
zero using the function m(T) oc (T c - T) f3 . For the DCA 
this extrapolation is shown by the solid line in the inset 
to Fig.|| for N c = 4. The values for T c obtained from the 
calculation in the symmetry broken phase and in the un- 
broken phase must agree, since as we have shown above, 
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FIG. 2: The transition temperature T c for the DCA (circles) 
and MCA (squares) when U — W = 1 versus the cluster 
size N c . For all values of N c the DCA prediction is closer to 
the exact result (T c — 0). Inset: Order parameter m(T) and 
inverse charge susceptibility x(Q) _1 versus temperature. The 
solid (dashed) line represents a fit to the functions m(T) <x 
(T c - T) with (3 = 0.245 (x(T) oc (T — T c ) -7 with 7 = 1.07). 



both the DCA and MCA are ^-derivable. This is illus- 
trated in Fig.| for the DCA. 

A comparison of the DCA and MCA estimate of T c is 
presented in Fig. |^. T c obtained from MCA (squares) 
is larger than T c obtained from DCA (circles). Moreover 
we find that the DCA result seems to scale to zero almost 
linearly in 1/N C (for large enough N c ), whereas the MCA 
does not show any scaling form and in fact seems to tend 
to a finite value for T c as N c — » 00. This striking differ- 
ence of the two methods can be attributed to the different 
boundary conditions. The open boundary conditions of 
the MCA cluster result in a large surface contribution so 
that Tmc > ^dc- This engenders pronounced mean field 
behavior that stabilizes the finite temperature transition 
for the cluster sizes treated here. For larger clusters we 
expect the bulk contribution to the MCA free energy to 
dominate so that T c should fall to zero. 

Complementary results are found in simulations of 



finite-sized systems. In general, systems with open 
boundary conditions are expected to have a surface con- 
tribution in the free energy of order 0(1/L c )[ff^. This 
term is absent in systems with periodic boundary con- 
ditions. As a result, simulations of finite-sized systems 
with periodic boundary conditions converge much more 
quickly than those with open boundary conditions Ju|. 

Summary. By defining appropriate Laue functions, 
we provide microscopic diagrammatic derivations of the 
MCA and DCA. We show that they are 4>-derivable, and 
that the lattice free energy is optimized by equating the 
irreducible quantities on the lattice to those on the clus- 
ter. The MCA maps the lattice to a cluster with open 
boundaries and consequently, the cluster violates trans- 
lational invariance. In contrast, the DCA cluster has 
periodic boundary conditions, and therefore preserves 
the translational invariance of the lattice. This differ- 
ence in the boundary conditions translates directly to 
different asymptotic behaviors for large clusters N c . As 
we find analytically as well as numerically, the surface 
contributions in the MCA lead to an average hybridiza- 
tion f of the cluster to the mean field that scales like 
1/L C as compared to the 1/L% scaling of the DCA. Since 
r acts as the small parameter for these approximation 
schemes, the DCA converges much more quickly than 
the MCA. These effects are more pronounced near a tran- 
sition, where the large surface contribution of the MCA 
stabilizes the mean-field character of the transition. Con- 
sequently, the DCA result for the transition temperature 
T c of the ID symmetric FKM model scales almost like 
1/N C to the exact result T c = 0, whereas the MCA result 
converges very slowly. Since the origin of this difference 
lies in the different boundary conditions we expect this 
primacy of the DCA over the MCA to hold generally for 
any model of electrons moving on a lattice. 
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